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Abstract. After a brief review of the production dynamics of prompt single photons in hadron col- 
lisions, I summarize a new QCD calculation of the transverse momentum distribution of continuum 
prompt photon pairs produced by QCD subprocesses, including all-orders soft-gluon resummation 
valid at next-to-next-to-leading logarithmic accuracy. Resummation is necessary to obtain reliable 
predictions, as well as good agreement with data from the Fermilab Tevatron, in the range of trans- 
verse momentum where the cross section is largest. Predictions are made for the Large Hadron 



Collider where the QCD diphoton continuum is shown to have a softer spectrum in transverse mo- 



mentum than the Higgs boson signal. 
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INTRODUCTION 



The study of prompt photons carrying large values of transverse momentum has a history 
almost as old as quantum chromodynamics (QCD). It ranks with deep-inelastic-lepton 
scattering, mas sive-lepton-pair production (the 'Drell-Yan process'), and jet production 
^ ■ as an important probe of short-distance hadron dynamics 1 . The appeal of prompt photons 

is that they are point-like, colorless probes of the dynamics of quarks and gluons, ones 
Qh| that escape unscathed through the colored medium of the high-energy collision. One 

motivation, in spin-dependent and spin-averaged scattering, is the access that prompt 
photons offer to the spin-dependent and spin-averaged gluon densities of hadrons. One 
of the two leading-order partonic direct subprocesses feeds directly from the gluon 
parton density through the gluon 'Compton' subprocess, qg — > yq. 

In reality, the promise and naive simplicity of prompt photons are compromised by 
two issues: fragmentation contributions and photon isolation. Fragmentation is a long- 
distance process in which a hard photon brems off a final-state quark (or gluon). The 
photon emerges as part of a jet if the opening angle between the quark and photon is too 
small. This 'showering' process must be parameterized by a non-perturbative single- 
photon fragmentation function Dy(z,fif) ■ Variable z is the fraction of the parent parton's 
momentum that the photon retains, and fragmentation scale \if serves as a boundary 
in momentum \x space between the long-distance dynamics included in Dy(z,fJ-f) and 
the short-distance next-to-leading order dynamics in the region /i > jif (i.e., large angle 
separation between the final photon and remnants of the fragmenting parton). Note that 
a new variable \Xf enters the mix along with the usual renormalization and factorization 



A photon is said to be prompt if it cannot be said to originate from the decay of a hadron, such as a 7T 1 
or 7], itself produced with large transverse momentum. 



scales \1r and /If- 

A theory calculation is most reliable for an inclusive observable since infra-red and 
collinear singularities can be shown to cancel between the real-emission and loop dia- 
grams in an inclusive calculation that is carried beyond the leading order in perturbation 
theory. However, Tevatron and Large Hadron Collider (LHC) experiments measure iso- 
lated photons in which part of the real-emission final-state phase space is excluded, 
resulting in incomplete cancellation of singularities in an analytic calculation yj]. Iso- 
lation is required since, otherwise, the prompt photon signal would be overwhelmed by 
secondary y's from hadron decays, and it is hard to discern a prompt y if it is buried 
in a jet (unless z is very large). In isolation procedures applied by the experiments, the 
energy of the hadronic remnants is required to be less than a threshold isolation energy 
Ej° in a cone AR = a/ At] 2 + A<p 2 around each photon, with At] and A<p being the sep- 
arations of the hadronic remnant(s) from the photon in the plane of pseudo-rapidity T] 
and azimuthal angle <p. The values of Ej° and AR characterize the measurement, and, in 
a theory calculation, the size of the fragmentation contributions depends on the assumed 
values of these parameters 2 . Theoretical representations of isolation are at best approx- 
imate. Calculations tend to require Monte Carlo modeling of fragmentation in order to 
implement experiment-like isolation selections. 

Calculations exist of prompt single photon production at next-to-leading order in 
QCD for both the direct and the fragmentation contributions. One approach is encoded 
in the Monte Carlo program JETPHOX An impressive comparison with experiment 
may be found in the DO paper [3]. 

Much more can be written about recent research on inclusive and isolated single 
photon production, but I turn to production of photon pairs in the remainder of this 
brief report, summarizing a new calculation of the diphoton cross section in perturbative 
QCD [4]. The motivation for this study comes, in part, from the fact that a Higgs boson 
with mass between 115 and 140 GeV may be identified at hadron colliders through its 
decay into a pair of energetic photons. Theoretical predictions of the QCD continuum 
'background' may be of substantial value in aiding search strategies. Moreover, the 
QCD calculation of continuum photon-pair production is of theoretical interest in its 
own right, and data from the Tevatron collider ||5J] offer an opportunity to test results 
against experiment. 

I focus on the transverse momentum (Qt) distribution of a pair of photons from 
the QCD continuum and from the Higgs boson signal. This distribution is important 
because its behavior affects the precision of the determination of the vertex from which 
the Higgs boson emerges. Greater activity in Qt associated with Higgs boson production 
allows a more precise determination of the event vertex. The expected shape of do / dQj 
can affect triggering and analysis strategies. Event modeling, kinematic acceptance, and 
efficiencies all depend on Qj. Finally, as is shown here, selections on Qt may be useful 
to enhance the signal to continuum ratio. 



In the case of production of photon pairs, discussed below, the two photons must also be separated in 
the r\ — (p plane by an amount exceeding the approximate angular size ARyy of one calorimeter cell. 



pp(qq,qg,gg) ~> YY-X. ^ = 1.96 TeV 
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FIGURE 1. Transverse momentum distributions of continuum diphotons for pp — > yyX at yfs — 
1 .96 TeV. The dotted histogram shows the result of the next-to-leading order calculation and the solid 
histogram the result after resummation at NNLL accuracy. The data are from the CDF measurement. 



CALCULATION OF THE DIPHOTON CONTINUUM 

Direct production of continuum photon pairs occurs from qq, q g, and gg scattering. In 
Ref. [4], next-to-leading order (NLO) perturbative cross sections are included, i.e., cross 
sections of €?(a s ) in the qq and qg channels, and G{a? s ) in the gg channel. Singular log- 
arithms arise in the NLO cross sections when Qj of the 77 pair is much smaller than its 
invariant mass Q. To describe properly the behavior of the Qj distribution in this region 
Qt < Q, where the cross section is greatest, all-orders resummation of initial-state soft 
and collinear logarithmic contributions is included, based on the Collins-Soper-Sterman 
(CSS) resummation procedure [6]. The logarithms are (re)summed into a Sudakov ex- 
ponent (composed of functions A(/i) and 5(/i)) and convolutions of the conventional 
parton densities f a (x,HF) with Wilson coefficient functions C. We include the ff(oc^) 
expressions for A(/i), &{ci^) expressions for B(fi), and the C-functions at order a s , for 
all subprocesses. These enhancements elevate the accuracy of the resummed prediction 
to the next-to-next-to-leading-logarithmic (NNLL) level. 



COMPARISON WITH TEVATRON DATA 

Our analysis provides the triple-differential cross section do / dQdQjd&ty (A0 is the 
difference of the angles of the two photons in the transverse plane). The calculation is 
especially pertinent for the transverse momentum distribution in the region Qj < Q, for 
fixed values of diphoton mass Q. It would be best to compare our multi-differential 
distribution with experiment, but the collider data tend to be presented in the form 
of single-differential distributions in Q, Qt, and A0, after integration over the other 
variables. For calculations appropriate at Tevatron energies, we impose cuts \yJ < 0.9 on 
the rapidity of each photon, and p T > p Tmin = 14 (13) GeV on the transverse momentum 



of the harder (softer) photon in each 77 pair. We choose Ej° = 1 GeV, AR = 0.4, and 
AR rr = 0.3. 

The Qj distribution is shown in Fig.[T] The resummed Qj distribution is well-behaved 
as Qj — > 0, unlike its fixed-order counterpart which is singular in this limit. The result 
after NNLL resummation is in good agreement with the absolute rate and overall shape 
of the data from the Collider Detector at Fermilab (CDF) collaboration at pp collision 
energy v/s= 1.96 TeVfl. 

In the two highest-<2r bins, there is some evidence of an excess experimental rate, 
perhaps a 'shoulder' in the data. Our calculations show that most of the shoulder events 
populate a limited volume of phase space characterized by small A(j) and Qj > Q. From 
a theoretical point of view, when Qj > Q, as in the shoulder region, the calculation 
must be organized in a different way ||7|] in order to resum contributions arising from the 
fragmentation of partons into a pair of photons with small invariant mass, a new feature 
of diphoton production with respect to single photon production. Adequate treatment of 
the light 77pairs is missing in our calculation and in all other calculations at present. This 
interesting region warrants further theoretical investigation. Our theoretical treatment is 
most reliable in the region Qj < Q. When this selection is made, the contributions from 
small A0 are efficiently suppressed, and dependence on tunable isolation parameters and 
factorization scales is reduced. After the selection Qj < Q, we expect that the large Qj 
shoulder will disappear in the experimental Qj distribution. 

Predicted distributions in the invariant mass Q of the photon pair and the separation 
may be found in Ref. [4]. Suggestions are made there for more differential analyses 
of the Tevatron data that would allow refined tests of the calculation. An important 
consequence of the resummation formalism is a logarithmic dependence on the diphoton 
invariant mass Q. We encourage the CDF and DO collaborations to verify this predicted 
broadening of the Qj distribution as Q increases. 

PREDICTIONS FOR THE LHC AND COMPARISON WITH THE 

HIGGS BOSON SIGNAL 

Several predictions are presented in Ref. 0] for the photon-pair continuum in pp col- 
lisions at s/s = 14 TeV. In this summary, I limit my discussion to differences expected 
between the spectrum of photon pairs from the Higgs boson signal and from the QCD 
continuum. The dominant production mechanism for the Higgs boson is gluon fusion, 
gg — > H — > 77, and we calculate the signal with the same order of precision in the QCD 
contributions as for the continuum. We include fixed-order initial state QCD corrections 
at 0(af) (NLO) and resummation at NNLL accuracy. These contributions are coded in 
the same numerical program used to compute the continuum, and we apply the same cuts 
(for example, on the values of pj and r\ of each photon) for both the signal and the con- 
tinuum. Use of one code minimizes potential deviations than can arise from codes with 
different physics content, parton densities, and numerical implementations. We compute 
the continuum in the range 128 < Q < 132 GeV and the signal at a fixed Higgs boson 
mass ma = 130 GeV. We use [Ir = 1±f = Q for the signal and the continuum. 

The cross section times branching ratio for the Higgs boson signal is substantially 
smaller than the continuum rate. We present distributions normalized to the respective 
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FIGURE 2. Comparison of the normalized transverse momentum distributions for the Higgs boson sig- 
nal and diphoton continuum at the LHC, both computed at NNLL accuracy. The continuum is calculated 
for 128 < Q < 132 GeV, and the Higgs boson mass is taken to be iuh = 130 GeV. 

total rates in Fig. El The signal and continuum peak at about 12 and 5 GeV, respectively. 
The dominant Sudakov exponent C> = = Aqg that controls gluon radiation for the 

qq and qg + qg initial states is less than Agg = Ca that controls the gg case. More 
gluon radiation in the gg case explains the broader Qj spectrum for the Higgs boson 
signal. The comparison in Fig. suggests that the signal to background ratio would be 
enhanced if a cut is made to restrict Qj > 10 GeV. Additional differences between the 
signal and continuum are shown in Ref. 0], including the rapidity difference of the two 
photons, the scattering angle in the Collins-Soper frame twah^M—yT) /2) = cos 0*, and 
the difference between the azimuthal angles of the photons. 
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